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ABSTRACT 

12 Bootis is a spectroscopic binary whose visual orbit has been resolved by interfer- 
ometry. Though the physical parameters of the system have been determined with 
an excellent precision, the theoretical modelling of the components is still uncertain. 
We study the capability of solar-like oscillations to distinguish between calibrated 
models of the system obtained by including in the stellar modelling different mixing 
processes. We consider different scenarios for the chemical transport processes: clas- 
sical overshooting, microscopic diffusion and turbulent mixing. For each of them we 
calibrate the stellar models of 12 Boo A and B by fitting the available observational 
constraints by means of a Levenberg-Marquardt minimization algorithm, and finally, 
we analyze the asteroseismic properties of different calibrated models. Several solutions 
with 12 Boo A in (or close to) post-main sequence and 12 Boo B on main sequence 
are found by assuming a thickness of the overshooting layer between 0.06 and 0.23 
the pressure scale height. Solutions with both components on the main sequence can 
be found only by assuming an overshoot larger in the primary than in the secondary, 
or a more efficient central mixing for 12 Boo A than for 12 Boo B. We show that the 
detection of solar-like oscillations expected in these stars would allow to distinguish 
between different scenarios and provide therefore an estimation of the overshooting 
parameters and of the properties of extra-mixing processes. 

Key words: stars: binaries - stars: oscillations ~ stars: interiors ~ stars: fundamental 
parameters - stars: individual: 12 Boo 



1 INTRODUCTION 

In standard stellar modelling convection is described by 
means of a local theo ry such as the Mixing Length Theory 
jBohm- VitenseHTosi '). and the extension of the convective 
regions is determined by the classical Schwarschild criterion 
(Vrad = Vad). This is One of the well known shortcomings of 
stellar structure and evolution theory. For intermediate and 
high mass stars, where a convective core is developed, the 
mass of this central mixed region plays a fundamental role in 
determining the lifetime and the luminosity of the main se- 
quence phase, and affects other relevant aspects such as the 
s-process nucleosynthesis o n the asymptotic giant branch 
l|Ventura fc D'AntonalboOSl '). the chemical evolution of the 
interstellar medium, and the size and composition of white 
dwarfs (see, e.g. [Chiosi 1998). 

We mean here by overshooting or extra-mixed region 
the mixing of material beyond the formal boundary of con- 
vection region set by the Schwarschild criterion. There is 
quite a large number of observational evidences indicating 
that this phenomenon of overshooting exists. For instance, 
com parisons between observations (binaries and open clus- 
ters, lAndersen et~al]|l990l : iRibas et al.ll2000l ) and theoreti- 



cal models indicate that stars have larger convective cores 
than predicted by theory. Moreover, numerical simulations 
and laboratory experiments suggest the overshoot of con- 
vective elements in the adjacent stable regions. A great 
effort has been done during the last 30 years to describe 
turbulent convection in stellar interiors modelling, how- 
ever, given the difficulty in estimating the rate of dissipa- 
tion of turbulent kinetic energy, it is still not possible to 
predict the thickness of the extra-mixed region from first 
principles. Therefore, the thickness of the overshooted re- 
gion (Aov) is usually parameterized in terms of the lo- 
cal pressure scale height Hp (Aov = aov Hp). The value 
of this parameter has been estimated by fitting theoret- 
ical isochrones to the observed color-magnitude diagrams 
of open clusters, and by modelling double-lined eclipsing 
binaries with we ll determined masses and radii (see e.g. 
iRibas et al1l2000l ). These calibrations also indicate that the 
exact value of the overshoot parameter might depend on 
the stellar mass: values of the order of 0-0.1 are required for 
masses between 1.1 and 1.5 Mq, and values fro m 0.2 to 0. 5 
are needed in the 1.6 to 9 Mp, interval (see.Ribas et al.ll2000l ). 
Moreover, as noticed in IVandenBerg fc StetsonI (|2004l ') and 
IVandenBerg et al] (|2006l ) . the overshooting parameter has 
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probably a steep dependence on mass in the range 1.3 < 
M/Mq < 1.55. 

Not only the exact amount of overshooting and its de- 
pendence on the stellar mass, but also the physical pro- 
cess responsi ble for extra-mixing, are still matter of a 
lively debate ("Macdcr fc Zahn'tgoj; lYoung fc Arnettll2005l : 
iPopielski fc Dz iembowski 2005). 

In this context the binary system 12 Bootis represents 
a valuable observational test of stella r modelling. Recent 
interferometric measurements allowed iBoden et al.l (|2005l ) 
to determine the masses of the components of this double- 
lined spectroscopic binary with a precision of about 0.3%. 
The masses of the components are very similar (mass ra- 
tio q ~ 0.97), but their luminosities are quite different. 
A preliminary determination of stellar parameters by corn - 
parison with theoretical isochrones (see lBoden et al.]|2005h . 
suggested that secondary component is still in the central- 
hydrogen burning phase, while the primary is more evolved 
and burning hydrogen in a shell. Given the lifetime of the 
sub-giant phase derive d from theoretical mo dels, this solu- 
tion is quite unlikelv ijAnderseii et al.|[l990l ). Furthermore, 
as shown bv lBoden et al.l "( 20051 ). the evolutionarv state of 
the fitted models is highly dependent on the model details 
(and in particular on the amount of core overshooting) and 
needs further investigation. The masses of 12 Boo compo- 
nents (1.416 and 1.372 Mq) are in fact in the transition re- 
gion, where the thickness of the "overshooting" layer seems 
to dep end on the stellar mass, metallicity and evolutionary 
stage (|VandenBerg et al.llioog '). 

Even if the masses are very precisely known, and both 
the components have the same age and initial chemical 
composition, there is a large number of possible parameter 
sets able to fit the "classical" observables. In this paper we 
present different scenarios based on different theoretical pre- 
scriptions for extra-mixing processes, which provide stellar 
parameters compatibles with observational constraints, but 
with different evolutionary state and internal structure. The 
classical observables do not allow to distinguish between dif- 
ferent solutions, and hence additional and independent ob- 
servational constraints are needed. In this paper we propose 
solar-like oscillations as additional observational constraints. 
Fortunately, the components of 12 Boo are located in a re- 
gion where solar-like oscillations arc expe cted to be exc ited 
(see e.g. iK icldscn & Bedding 1995. Houdek et al.lll999l and 



ISamadi et al.ll2005l ). 

Solar-like oscillations have been ob served in several sin- 
gle s ta rs (e.g. rj Boo t is, an d fj, Arae, iKieldsen fc Beddind 
1 19951 : iBouchv et al.l |2005|. re s pectiv ely) and binary 
systems (Procyon, iMartic et al.l |2004| ) as well as in 
both components of th e visu a l binary system a Cen- 
tauri ([Bouchy fc C arrier '2002": Carrier fc iBourbanI l2003l : 
2004 : Kjeldson ct al. 2005). In the latter, as 



tauri l l ijoucny 
iBedding et al.l 



shown by iMiglio fc MontalbanI (|2005l ) , the combination of 
asteroseismic and precise "classical" constraints (masses, 
radii, luminosities . . . ) significantly improves the determi- 
nation of the system fundamental parameters. While the 
study of a Centauri provides a test of our knowledge of stel- 
lar structure in conditions that are slightly different from the 
Sun, a detailed investigation of 12 Bootis will be a relevant 
step in understanding the structure of the central regions 
for this particular area in the Hertzprung-Russell diagram 



(HRD) where the "transition" from radiative to convective 
core models takes place. 

In Sect. [2] we summarize the observational knowledge 
about 12 Boo, and we describe the constraints will be used 
in the caUbration process. The latter is described in Sec. |31 
Among the free parameters, the most important is the exten- 
sion of the overshooting region. We will also analyse possible 
solutions by taking into account different physical processes 
able to change the thickness and features of the extra mixed 
region. In Sec. Uwe study how the differences in the stellar 
structure for the models fitted in different scenarios are re- 
flected on the oscillation frequencies, and we suggest which 
kind of observational data could give us enough information 
to constraint the evolutionary state of 12 Boo, and, possi- 
bly, also the properties of the transport processes (if any) in 
the stellar central region. Finally, in Sec. [S] we present our 
conclusions. 



2 12 BOOTIS: OBSERVATIONAL 
CONSTRAINTS 

With the rapid development of optical interferometry in the 
last years, closer and closer binary systems have become 
resolvable. Thanks to interferometric orbits an increasing 
sample of double- lined spectroscopic binaries (SB2), includ- 
ing non-eclipsing systems, can be used for accurate stellar 
parameter determination. Interferometry essentially yields 
the orbital inclination that, combined with accurate radial 
velocity measurements, can provide mass determination at a 
level better t han 1%, i.e. suita ble for strict stellar modelling 
tests (see e.g. lAndersen|[l99]] l . 

A good example of an ideal target for this technique 
is 12 Boo, a bright (V=4.83) non-eclipsing SB2 with an or- 
bital period of about ten da ys (P = 9.6) and a (composite) 
spectral type F8 IV -F9 I VwiBarrvl (Il970l'l (w here "w" stays 
for weak metallic lines). iBoden et al. (|2005h obtained long 
baseline interferometry of 12 Boo with the Palomar Testbed 
Interferometer (PTI) and the Navy Prototype Optical In- 
terferometer (NPOI). They secured, as well, new high res- 
olution (ec helle) spectra, as t h e previous spectroscopic ob - 
servations (|Abt fc Levvlll976l : [Pe Medeiros fc Udrvlll999l ). 
were the major limit to the precision of derived parameters. 
The interferometric visibility data and the radial velocity 
curves were simultaneously solved to derive the orbital pa- 
rameters. The components could not be resolved (the ap- 
parent orbital semiaxis is only 3.45 mas), so that the indi- 
vidual radii had to be derived by indirect means, namely 
by t he Infrared Flux Method by Blackwell and collabora - 
tors (iBlackwell et al ] 1 19901 : IBlackwell fc Lvnas-Gravlll994l ). 
The information about the total and the individual bolo- 
metric fluxes (needed for radius determination) was derived 
by a fit of multi-wavelength archival data and from the in- 
band intensity ratio of the interferometric measurements. 
Furthermore, the analysis of the observed spectra provided 
an estimation of the component effective temperatures and 
of the luminosity ratio. A further improv ement in the de- 
termin ation of the spectroscopic orbit by iTomkin fc Fekell 
(|2006l ). combined with the same interfe rometric data, con - 
firms the values of the masses derived bv lBoden et~al] (|2005D 
and reduces the associated error bars. The relevant values 
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Table 1. Orbital parameters and magnitude differences for 12 
Boo 



Table 2. Observational constraints as adopted in the modelling 



Parameter 


value 


reference 


P(d) 


9.6045529 ±4.8 IQ-'^ 


1 


e 


0.19268 ± 0.00042 


1 


i 


107?99 ± 0.077 


2 


a (mas) 


3.451 ± 0.018 


2 


AiCciT (mag) 


0.589 ± 0.005 


2 


A/fciT (mag) 


0.560 ± 0.020 


2 


AV (mag) 


0.560 ± 0.020 


2 


Ka (Km s-i) 


67.286 ± 0.037 


1 


Kb (Km s"!) 


69.30 ± 0.050 


1 



" References: It lTomkin fc Fekell l|2006l ).2') iBoden et al.l l|2005h 



and their uncertainties are reported, for the reader's sake, 
in Table [U 

Other relevant information on the system are its near- 
solar chemical composition, with a mean value of [Fe/H] = 
0.01 ± 0.08, derived from Stromgren photom etry and de- 
tailed spectroscopic ana l ysis (Duncan 1981: Balachandranl 
Il990l : iLebre et al.l Il999l : iNordstrom et al.1 |2004| ) and the 
component rotation al velocities of 14.0 and 12.0 km s~^ 
jBoden et al.1 |2005| ). The latter are slightly faster than 
the expected corotation rate of, respectively, 12.4 and 
9.3 km s-^ 

The resulting global picture of 12 Boo is that of a bi- 
nary with very similar components (masses, effective tem- 
peratures, rotation) except an unexpected magnitude differ- 
ence: both spectroscopic and interferometric measurements 
indicate a difference of 0.48-0.58 mag, depending on the 
band. That is quite hard to explain in a system with co- 
eval and non-interacting components (both of them are well 
inside the critical Roche lobe). 

The high precision of the radial vel ocity measurements 
(the typical error is 0.5 km s~^ in iBoden et al. I l2005l : 
iTomkin fc Fekelll2006l ) and the short orbital period rules out 
the possibility of making the masses unequal by "hiding" a 
small fraction of a component mass in a third unseen object. 
Unless rather improbable configuration of the third object 
orbit are assumed (i.e. orbit almost normal to the line of 
sight), even an object of a few hundredth solar masses or- 
biting one component would easily be detected from the 
reflex motion of the more massive companion. Besides, such 
a planetary orbit can be excluded on the ground of stability 
consideration, see e.g iHolman fc Wiegertl l| 19991 ) and refer- 
ences therein. 

The observational constraint s we assume i n our mod- 
elling are taken from Table 5 of iBoden et al. I l|2005l ) that, 
for the sake of convenience, we repeat here in Table [H As 
the above mentioned spectroscopic and photometric studies 
of 12 Bootis indicate a metallicity within 0.1 dex of the so- 
lar value, given this uncertainty and a further uncertainty 
on the value of the solar metall icity, we adopted the solar 
Z/X of [g revesse fc Noels! (119931 ) considering a conservative 
error bar, as reported in Tabled 





A 


B 


M/M0 


1.4160±0.0049 


1.3740±0.0045 


Teff (K) 


6130±100 


6230±150 


L/Lq 


7.76±0.35 


4.69±0.74 


Z/X 


0.0245±0.01 


0.0245±0.01 



3 MODELLING 12 BOOTIS 

In the following we will use the term "calibration" to de- 
note the process of determining the stellar model param- 
eters, such as initial composition and age (both the same 
for A and B components), satisfying - within 1-a - the ob- 
servational constraints: Ma, Mb, Z/X, La, Lb, Teff. a and 

Tcff.B- 

As both system components are massive enough to de- 
velop a convective core during the main-sequence and no a- 
priori assumption on their evolutionary state can be made, 
the dependence of T^s and luminosity on the stellar param- 
eters (e.g. a g e and overshooting) is highly non-linear. As 
iBoden et al.1 l|2005l ) already pointed out, the result of sys- 
tem fitting strongly depends on the model details, and, in 
particular, on the amount of convective overshooting in the 
core. 

The aim of this paper is to study the possible solu- 
tions (not just to find a single solution) and to verify the 
ability of discriminating among different scenarios. We ap- 
ply, theref ore, the minimization technique adopted and de- 
scribed in iMigho fc MontalbanI (|2005l ) with the purpose of 
determining the chemical composition and age for each value 
of the overshooting parameter (assumed to vary from 0.0 to 

0. 50 in steps of 0.01). The parameters left free are: the ini- 
tial chemical composition (described by the hydrogen and 
heavy-elements mass fraction, X and Z respectively), the age 
and the component masses. These are adjusted in order to 
fulfill the observational constraints in Table[2] A standard 

is used as goodness-of-fit estimate, so that each observable 
is taken into account with its uncertainty. We did not allow 
omlt to be adjusted and fbced its value to qmlt = 1.84, 

1. e. the value provided by the Solar calibration. This choice 
is justified by the fact that both stars are quite close in the 
HRD and that different values of OfMLT are not expected. Be- 
sides, a variation of Qmlt implies a change of Tch that can 
be easily balanced by a change of the chemical composition. 

The search for a solution fulfilling the observational con- 
straints on the HR diagram is usually biased by an addi- 
tional criterion, though this is not included in the x^'- ^ 
preference towards long lasting evolutionary phases, as the 
probability of finding a star in a given evolutionary phase 
depends on how "rapid" that phase is. Evidently, it is less 
likely to observe a star in its second gravitational contrac- 
tion than on the main sequence. This is however a statis- 
tical argument, that can be verified on large samples of 
stars (e.g. in clusters), but cannot be applied to a single 
case. In addition to that, the duration of each evolutionary 
stage depends on whether overshooting is considered i n the 
models or not. A s pointed out e.g. by iMaedeij (Il975l ) and 
iNoels et al.l l|2004h , the duration of the thick-shell- hydrogen 
burning phase depends on the amount of overshooting as 
the latter increases the mass of the isothermal helium core 
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built during the core hydrogen-burning phase. As an exam- 
ple, a I.4M0 star without overshooting s pends in tha t phase 
~ 20% of the main-sequence lifetime (see llbenlflQQll . for in- 
stance). An overshooting parameter aov — 0.1 implies a 
decrease of the duration of this evolutionary phase down to 
4% of the main-sequence lifetime, i.e. lasting as long as the 
the overall contraction phase. In our modelling we do not 
make, therefore, any a-priori assumption on the particular 
evolutionary phase of the solution. 



Stellar Models 

The stellar model sequences are computed with the code 
CLES (Code Liegeois d'Evolutio n Stellaire). The op acity ta- 
bles are those of OPAL96 (I glesias fc Rogerslll996l) comple - 
mented at T < 10000 K with I Alexander fc FergusonI (|l994l ') 
opacities. The metal mixtur e used in the opac i ty tab les, is 
the solar one according to iGrevesse fc Noeld (I1993I ). The 
nuclear ene rgy generation routi n es are based on the cross 
sections by ICaughlan fc Fowler! l| 19881 ) and updated with 
the recent measuremen t of the ^'^ N{p,^y^O reaction rate 
jFor micola et al."2004''). The weak screening factors come 
from Salpctcr (1954), and the equation of state used i n 
the computations is OPALOl (|Rogers fc Navfonovl |2002| ). 
Convection trans port is treated with the classical mixing 



Table 3. Parameters of the models. 



len gth theory b y Bohm-Vitensel l|l958l ) in the formulation 
bv lCox fc Giulil (11 968') ; and atmospheric boundary condi- 
tions given bv lKurucz (,1998. ) are applied at T = T^s- 

Chemical mixing in the formal convective region and 
in the overshooting layer is treated as instantaneous. Be- 
sides, the thickness of the overshooting layer, Aov, is ex- 
pressed in terms of a parameter aov, so that Aov = 
Qov X min(rcv, Hp{rcc)), where rcc is the radius of the con- 
vective core. The temperature stratification in this over- 
shooting layer could be radiative or adiabatic (e.g. Zahn 
1991). At any rate the effect of this thin layer on the evolu- 
tionary track is very small, in our models a radiative strati- 
fication has been assumed. 

We have computed as well models that, in contrast with 
the assumed instantaneous chemical mixing, include a dif- 
fusive treatment of mixing in the formal convective region 
and in the extra-mixed l ayer. These computatio ns were done 
with the ATON3.0 code P'Antona et al.ll2005l ) for the same 
physics: MLT, OPAL1996 and Alexander fc Ferguson opac- 
ities, as well as OPAL equation of state. 

We have also studied (see Sect. I3.2"|l the effect of other 
chemical mixing processes, such as diffusion and gravita- 
tional settling of helium and heavy elements. Microscopic 
diffusion has b een implemented in CLES code following the 
formulation bv lThoul et"aLl (Il994l ). 

Finally, the effect of an additional turbulent mixing, 
such as the one expected to be induced by rotation, is ad- 
dressed in Sec. 13.31 





Qov 


^0 


Zo 


Age (Gyr) 


a 


0.06 


0.268 


0.0195 


3.08 


b 


0.15 


0.267 


0.0198 


3.31 


c 


0.23 


0.267 


0.0200 


3.58 


d 


0.37, 0.15 


0.281 


0.0186 


3.24 



Fig. [T]a: the primary component is close to the maximum 
luminosity in the shell-hydrogen burning phase, whereas 
the secondary, with a Xc=0.03, is at the end of its cen- 
tral hydrogen-burning phase. As aov increases, the evolu- 
tionary state of component A approaches the TAMS. With 
Qov — 0.15 (Fig. [T]b) the primary is burning hydrogen in a 
thick shell while component B is well on the main sequence. 
A further increase of the overshooting parameter aov ^ 0.2 
places the primary in the rapid overall contraction phase 
(Fig. [T]c). Still larger values of aov make Lb increase out 
of 1-CT and therefore worsen the fit. 

The observed luminosity ratio, La/Lb 1.65 (from a 
minimum of 1.35 to a maximum of 2.) is the main obstacle 
for a solution with both components on the Main-Sequence. 
In fact, given that both components have the same age and 
initial chemical composition, if the same kind of mixing 
processes are considered in the center, the luminosity ra- 
tio in the main-sequence is determined by the mass ratio 
{q ~ 0.97), and does not significantly depend on the choice 
of the other common model parameters {X, Z, aov). 

Regardless of the set (X, Z, aov) chosen, all the mod- 
els computed with aov. a = aov.B provided {La/ Lb)ms — 
1.15. A higher value of the luminosity ratio can be reached 
if aov, A > aov,B. Therefore, the only way we found to 
place both components in their MS phase is to assume a 
different efficiency of mixing processes in each component. 
In Fig. [T]d we show a calibration solution obtained by as- 
suming aov. A — 0.37 and aov,B — 0.15. In this case both 
12 Boo A and B are well in the MS. 

The values of stellar parameters collected in Table [3] 
are those used in the figures and comparisons in this pa- 
per. It has to be noticed, however, that the uncertainty on 
Z/X reflects on an undertainty of 0.02 in the initial helium 
mass fraction (Yq) and of 0.2 Gyr in the age of the models. 
Besides, the calibrations with aov, a ~ aov,B and similar 
initial chemical composition provide an age difference of the 
order of 15%, and the calibration with aov, a > aov,B im- 
plies an age in the range bracketed by the aov, a = aov,B 
models. Unfortunately, the uncertainty in Z/X is too large 
to constrain the models. 

All the different scenarios presented here cannot be 
clearly discriminated by means of the available observa- 
tional constraints, nevertheless, the internal structure of 
these models are significantly different, and these differences 
shall show up in their asteroseismic properties (see. Sec.|4]). 



3.1 Overshooting 

We performed several calibrations for different values of 
the overshooting parameter and we found solutions with 
aov ^ 0.06 (see. Table [5}. By changing aov the evolu- 
tionary state of the system changes. In detail: for the lowest 
values of overshooting the situation is that represented in 



Diffusive overshooting 

The classical treatment of overshooting is an extension of the 
instantaneously mixed convective core by a distance Aov. In 
the diffusive approach of convection, chemical mixing is de- 
scribed by the diffusion equation with diffusion coefficient 
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Figure 1. HR diagram showing solutions found with different values of the overshooting parameter: a) oqv = 0.06, b) aov = 0.15, c) 
aov = 0.23, and d) oqv A = 0.37 and aov B = 0.15. Model parameters are described in Table |3] Error boxes corresponding to 1- and 
2-(T error bars in L and T^g are represented by continuous and dotted lines. 



D = l/3vconvA, where Vconv is the average turbulent veloc- 
ity and A the convective scale length. In a MLT treatment 
Vconv oc (V — Vad) and A (X Hp, which imply t^conv = at 
the formal boundary of the convective region (V — Vad)- 
In order to define a continuous diffusion coefficient, also 
in the overshooting region, ATON code extrapolates the 
log?;conv vs. logP relation, to get the turbulent velocity at 
the boundary, starting from a point inside the convective re- 
gion where the pressur e is 5% larger than at the boundary. 
Furthermore, following IXiond |l985), turbulent velocity is 
assumed to exponentially vanish outside the convective re- 
gion, and the thickness of the mixed region is determined by 
the parameter f3 which describes t he exponential vanish ing 
of turbulent velocity. As shown by IVentura et all (|j^99^) , a 
/3 value ten times smaller than the aov in classical instanta- 
neous mixing produces the same effect on the HR diagram. 

With this formulation of overshooting we find the same 
set of solutions as with instantaneous mixing, this time for 
a parameter P in the range 0.004-0.035. This is of course re- 
assuring as a completely independent stellar evolution code, 
and a different description of overshooting have been used. 
The chemical composition profiles at the border of the con- 
vective core are, however, much smoother with the diffusive 
overshooting than with the classical treatment. 



3.2 Effects of diffusion 

iRichard et aD l|200lh . iMichaud etahl (|2004l ) and iRichardI 



i2005[) showed that microscopic diffusion can affect the size 
and evolution of a convective core. In those papers, thermal 
diffusion, gravitational settling and radiative accelerations 
were included in the models to derive the evolution of chem- 
ical distribution. CLES code does not take radiative acceler- 
ations into account. Nevertheless, there is no observational 
evidence in 12 Boo of surface chemical peculiarities that 
could be produced by gravitational-settling and radiative- 
acceleration, and the latter is not expected to play a mayor 
role in the central region (see e.g. Richard ct al. 2001). 

We do not expect a large effect of radiative acceleration 
in the 12 Boo mass domain, since the thickness of their con- 
vective envelopes is quite large. Even if at the beginning of 
the MS component A has a shallow convective envelope, and 
therefore the role of radiative accelerations could be non- 
negligible, as the star evolves its convective region becomes 
deeper. As a consequence, the effect of radiative acceleration 
is reduced, and the chemical composition of the convective 
region re-homogenized. We decided, therefore, to consider 
and compare models accounting for diffusion of H and He, 
and of H, He and Z. 

iMichaud et al.l l|2004 ) showed that microscopic diffusion 
can induce an increase of the convective core for a narrow 
range of masses, from 1.1 to 1.5 M0, and that the effect 
decreases rapidly with increasing stellar mass. We also note 
that in this mass domain, the convective core mass increases 
during the MS evolution instead of decreasing, as it occurs 
for larger masses. As a consequence, a sharp gradient of 
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10^ 2x10** SxlO** 

Age (yr) 



Figure 2. Fractional mass of tlie convective core as a function 
of the age in no diffusion models: without overshooting (solid 
line) and with Oov =0.1 (long-dashed line); and in no overshoot- 
ing models including microscopic diffusion of He and If (dashed- 
dotted line) and of He,H and Z (dashed line). 
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Figure 3. Evolutionary tracks of the models described in Fig. [2] 
in a log Tefi-Age diagram. 



chemical composition appears at the convective core border, 
making the diflFusion process much more efficient in that 
region. 

As shown in Fig. [5] for a 1.4 M© star, the increase of 
the core mass fraction when diffusion of H, He and Z is con- 
sidered is of only ~ 5%, and the age of the turnoff is only 
slightly inc reased (Fig. [51) . This is in agreement wit h the 
findings of iMichaud etaP ()2004l ) and iRichardI (j2005l ) who 
consistently included radiative accelerations in their com- 
putations while considering Z-diffusion. An obvious conse- 
quence is that, in presence of diffusion, the amount of over- 
shooting needed to find 12 Boo A and B in a given evolu- 
tionary state is reduced with respect to the values obtained 
in the previous section. In fact, the minimum amount of 
overshooting required is slightly smaller, i.e. aov — 0.04 
instead of 0.06. The evolutionary state of the calibrated sys- 
tem changes, with increasing aov, in a way similar to that 
described Sec. 13.11 

On the other hand, the increase of the opacity in the en- 
velope due to the settling of He and Z provides evolutionary 
tracks cooler than those without diffusion (this effect is even 
more important if only H and He diffusion is considered, see 
Fig. [3} . At the end of the main sequence the surface ZjX 
is ~ 10% (diffusion of H and He) or ~ 25% (diffusion of H, 
He and Z) smaller than the initial value, nevertherless, given 
the uncertainties in the observed Z/X, this implies only a 
slight change in the Z and X initial values in the calibrated 
models. 

In addition to a slight enlargement of the convective 
core, diffusion leads to a smoother chemical composition pro- 
file at the edge of the convective core (see Fig. llUp . Whether 
this property affects the oscillation modes of the components 
will be addressed in see Sec. |4l 



3.3 Other sources of extra-mixing in the core 

The shear instability induced by stellar rotation is consid- 
ered as one of the first causes of chemical mixing in the ra- 
diative interior of the stars. Different approaches have been 
proposed to treat the effects of rotati on on transport of 
angu l ar momentum and chemicals (see iMaeder fc Mevneg 
I2OOOI : iMathis et al] |2004| . for a review). Since our stellar 
evolution code does not include a consistent treatment of 
rotational effects, we simply include the chemical turbulent 
mixing by adding a turbulent diffusion coefficient (Dt) in 
the diffusion equation. In our parametric approach Dt is 
assumed to be constant inside the star and independent of 
age. These simple prescriptions are inspired by the results 
obtained bv lMathis et al.l (|2004l ) from a consistent study of 
rotational effects in a 1.5 M© star (see their Fig. 2). 

We have then calibrated our models with different val- 
ues of Dt. We find, similarly to the case of including over- 
shooting (see Sec. 13. ip . that we are able to fit the system 
at different evolutionary stages depending on the value of 
Dt; the values needed to fit the system are in the range 
30 — 170 cm^ s~^. We note that these values are of the same 
order of the total diffusion coefficient resulting from the cal- 
culat ions for a 1.5 M0 model as presented in iMathis et al.l 
(|2004l ). As already seen for overshooting, in order to fit the 
system with 12 Boo A and B in their main-sequence, a differ- 
ent efficiency of mixing must be assumed in each component: 
Dt,a = 330 cm^ s"\ and Dt.b = 120 cm^ s"^ 

The simplified parametric treatment of rotationally in- 
duced mixing used in this work has the aim of showing that 
if an extra-mixing process, different from overshooting, is 
acting near the core it will produce a different chemical com- 
position profile in the central regions of the star. Whether 
these different profiles are reflected or not on the solar-like 
oscillations features will be analysed in Sec l4.2l 
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4 ASTEROSEISMIC PROPERTIES 

Stellar models predict the presence of a convective envelope 
for the range of masses of 12 Boo A and B and, therefore, 
solar-like oscillations are expected to be excited in both com- 
po nents of the system. By m e ans of the scaling re l ations 
of iKieldsen fc Beddind (|l995l ), iKieldsen fc Beddind (1200 il l 



Houdek et 



and the theoretical predictions of 
and using the precisely determined values of the mass 
and the luminosity, we estimate that the modes with the 
largest amplitudes should be around fc'max ~ 700 /xHz and 
I'max =1200 /iHz in 12 Boo A and B, respectively. The peak 
amplitude in the power-spectrum is expected to be 3-4.5 
times larger than the solar-one for the A component, and 
2.5-3.5 times larger then the solar one for the B component. 

The components of 12 Bootis are massive enough to 
develop a convective core during main sequence. The com- 
bined action of nuclear burning and convective mixing is 
responsible for the development of a steep chemical com- 
position gradient at the boundary of the convective core. 
As a star leaves the main-sequence, the increasing central 
condensation, together with the steep chemical composition 
gradient in the central regions, leads to a large increase of the 
buoyancy frequency (A'') and, therefore, of the frequencies of 
gravity modes. The latter interact with pressure modes and 
affect the properties of non-radial solar-like oscillati ons by 
the so-called avoided crossin g phenomenon (see e.g. lOsakil 
1 19751 : lAizenman et al.l Il977t ) . The frequency of the modes 
undergoing an avoided crossing (also called mixed-modes) 
are therefore sensitive probes of the core structure of the 
star. 

As a first example let us consider the evolution of the 
frequencies of the model of component A corresponding the 
scenario b of Table |3] (hereafter A6). In Fig. |4]we plot the 
evolution of frequencies for £ = (dotted lines) and £ = 1 
(solid lines) modes, and show how y-modes, whose frequency 
strongly increases after the overall contraction (~ 3.27 Gyr), 
interact with pressure modes of same degree. This interac- 
tion results in non-radial modes with a mixed p— and g— 
character and frequencies that significantly deviate from a 
regular spacing between consecutive overtones. 

Since the thickness of the mode evanescent region 
(where N < ui < Se, with Se being the so-called Lamb 
frequency) increases with the degree i (see propagation di- 
agram in Fig. [SJ, the interaction between p and g modes 
decreases. As a consequence, in modes of degree £ > 1, the 
behaviour of gravity a nd pressure mode is generally bet- 
ter s eparated Jsee e.g. IChristensen-Dalsgaard et al.l 1 19951 : 
iMore l ct al. 200lh : g-modes tend to remain concentrated in 
the core and would hardly be detectable due to their sig- 
nificantly larger inertia. This is in fact the case for the our 
models A a and Ab where only £ — 1 modes present a mixed 
character. Nevertheless, for the less evolved models (Ac for 
instance) also £ — 2 mixed modes appear in the solar-like fre- 
quency domain. In the following discussion we will address 
mainly the properties of mixed modes. 

Given the MS status of 12 Boo B, we will first discuss 
the seismic properties of component A, and only in Sec. 
14.1.11 consider the additional information that the observa- 
tion of solar-like oscillations of component B would add to 
the modelling. 




Figure 4. Evolution in time of the frequencies of model Ab. Dot- 
ted and solid lines represent radial and £ = 1 modes with radial 
order n from 11 to 23. The frequencies of model Ab are high- 
lighted by a thick vertical line. Vertical dashed lines represent 
±1 (T interval in the observed Toff. 



100 




Figure 5. Propagation diagram for model Ab (see Table [S] 
Fig-CJ- The solid line represents the Brunt- Vaisala frequency A'^, 
the short- and long-dashed line the Lamb frequency Si for, re- 
spectively £ = 1 and £ = 2. The domain of solar-like oscillations 
falls within the horizontal dotted lines. 



4.1 Probing the evolutionary status 

As already investigat ed in the case of rj Bootis (see e.g. 
iDi Mauro et al.l |2004| ) , the frequencies of modes undergo- 
ing an avoided- crossing are direct probes of the evolutionary 
status of intermediate mass stars. 
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Figure 6. Large frequency separation as a function of the fre- 
quency for the models described in Fig. [T] Solid and dashed lines 
connect, respectively, Au for radial and £ = 1 modes. 

In Fig. [6] we compare the large frequency separatiorQ 
Av for radial and i = 1 modes corresponding to the models 
of component A presented in Fig. [T] Though radial modes do 
not give information on the evolutionary status, £ = 1 modes 
allow a clear discrimination among the scenarios. Modes of 
mixed p and g character are not present in the model on the 
main sequence (Fig. [Sid) and is almost constant, indicat- 
ing therefore a regular spacing between modes of consecutive 
radial order (both radial and £ = 1). On the other hand, as 
more evolved models (Aa, Ab and Ac) are considered, the 
avoided crossings effects become more evident. Therefore, for 
the model in the second overall contraction (Ac, Fig.|S]c) the 
interaction between p and g modes changes by up to 10 /xHz 
the modes with frequencies lower than 600 /xHz (correspond- 
ing order-mode, n, lower than 12). That can also be seen in 
Fig. [4] for models at 3.265 Gyr and the lowest n modes. For 
more evolved models (Aa and A6) the avoided crossing phe- 
nomenon clearly breaks the regular frequency spacing in the 
domain of expected solar-like oscillation (Figs.|6]a and |6lb). 

The radical variation of the properties of £ = 1 modes 
derives from a considerable change of the Brunt- Vaisala fre- 
quency near the center of the star (see Fig. [7)l as hydrogen 
is exhausted throughout the convective core and the star 
undergoes an overall contraction. As models of component 
A at different evolutionary stages are obtained for differ- 
ent amounts of core overshooting, the detection of mixed 
modes and their frequencies would allow to constraint the 
value of the overshooting parameter or the extension of 
the mixed central region. While model Aa (computed with 
ctov ~ 0.06) presents avoided crossings at v — 780 /xHz 
(n = 15) and u = 1100 /iHz (n =24) model Ab (qqv = 0.15) 

^ Computed as the difference between modes of same degree and 
consecutive frequencies. 



Figure 7. Propagation diagram in the core of the models for 
which large frequency separation has been plotted in Fig. |6] Solid 
line represents the Brunt- Vaisala frequency, and the dashed line 
the Lamb frequency for £ = I. 



does it at = 700 /iHz (n=13) and u = 950 /xHz (n=20). 
These differences between the frequencies of the modes pre- 
senting mixed p-g character (~ 100 piHz), as well as the 
differences of the frequency spacing between two consec- 
utive avoided crossings (~ 20/iHz), are large enough to 
be detected with the current observational techniques. Fur- 
ther more, 12 Boo A has a rotational velocity 14 km s~^ 
(Bo den et aT]|2005l ') that, assuming an inclination of the ro- 
tational axis equal to the orbital one (108°), implies a rota- 
tional splitting between {£ — l,m — —1) and (^ = l,m = 1) 
modes of only 2.5 /^Hz. 

Distinguishing between scenarios Ac (aov = 0.23) and 
Ad (qov = 0.37 for 12 Boo A and aov = 0.15 for 12 Boo B) 
on the base of £ = 1 mixed modes might be less straight- 
forward since the effect of avoided crossing occurs at lowest 
frequencies, where the amplitudes of the modes are predicted 
to be small. Therefore, no detection of avoided crossing could 
imply, either that 12 Boo A is well in MS, and has therefore 
a large central mixed region, or that it is close to the TAMS. 
Luckily, for models in the second overall contraction (SOC), 
the changes in stellar structure are such that £ = 2 mixed 
modes appear in the predicted solar-like spectrum, while for 
MS models (Ad) they would have very low frequencies. 

4-1.1 Seismology of component B 

Though in all calibrations we considered component B as a 
main-sequence object and, therefore, avoided crossings are 
not expected, the detection of solar-like oscillations in the 
secondary would, nonetheless, add constraints to the mod- 
elling of the binary system. A first additional constraint 
comes from the average value of the large frequency separa- 
tion {{Au)) that, for main-sequence stars, is well known to 
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Figure 8. Position in the HR, diagram of 12 Boo A and B with 1 
and 2-cr error boxes in logL and logT^ff. Also shown are linos of 
constant radius deduced from the (known) masses of the compo- 
nents and the large frequency separation. The large separation is 
assumed to be known with an uncertainty of 1 (continuous lines) 
and 2 /iHz (dotted lines). 
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Figure 9. The small frequency separation 5uq2 as a function of 
frequency in models of component B corresponding to the scenar- 
ios a — d described in Fig. [T] 



represent an estimate of the mean density and, hence, of the 
steUar radius, if the mass is known. The availab le estimate 
of component radii in 12 Boo (|Boden et al.ll2005l ) is based on 
the Infrared Flux method and yields: 7?a = 2.474±0.096 R© 
and 7?B = 1.86±0.15 R©. However, a more precise and inde- 
pendent determination of the stellar radii could be obtained, 
thanks to the small uncertainty on the masses (0.3%), from 
the knowledge of (A//). In particular for component B, the 
measurement of (A;/), would result in a significantly smaller 
uncertainty on its location in the HR diagram and, there- 
fore, in tighter constraints on the modelling. In Fig. [8] we 
show the uncertainty on the radii deriving from (A;/) values 
supposed to be known with an accuracy of 1 and 2 /iHz. 

Information on the central structure of component B 
could be provided by the small frequency separations Svq2- 
In Fig. [5] we show 5uo2 as a function of the frequency, as ex- 
pected for the models of 12 Boo B corresponding to the four 
scenarios of Fig. [T] We notice that the dependence of 5po2 
on u is mainly determined by the location of the disconti- 
nuity in sound speed derivative at the border of convective 
core, while the age of the model changes the mean value of 
5vq2- The spread of &V02 for the four considered scenarios 
increases with the frequency and is of the order of 2 /iHz at 
1500 /iHz and of 1.5 /iHz at Vmax- The age difference be- 
tween scenarios 6 and d (note that both have qqv = 0.15) 
implies a shift of only 0.5/iHz in 5uo2- 

Can 5vo2 of component B help to discriminate between 
scenarios c and d in absence of avoided crossing detection? 
Fig. |9] shows that the difference of 5vo2 at i^max is only 
0.2 /iHz and reaches 0.5 /iHz at 1500 /xHz. Furthermore, 
the frequency of ^ = 2 modes will be affected by rotational 
splitting (expected to be ~ 1 — 1.5 fJ^Hz), so that an accurate 
and reliable measurement of the small separation (expected 
to be around 3 /xHz) will be a rather difficult task. Other 
frequency combinations, such as the small separation doi are 
known to be very sensitive to the core structure of main se- 



quence stars (see e.g. [Roxburgh fc Vorontsovl2003l and refer- 
ences therein), but the larger number of frequencies needed 
to compute these small differences, increases even more their 
observational uncertainty. 

4.2 Probing the chemical composition gradient in 
the core 

Is there any other information that the knowledge of the fre- 
quency of an avoided crossing can give us? The appearance 
of avoided crossings is due to the evolution of frequencies of 
gravity modes and their interaction with acoustic modes of 
similar frequency. The frequency spectrum of gravity modes, 
on its turn, is mainly determined by the behaviour of A'' in 
the central regions of a star {u (x J N/rdr). Furthermore, 
the mean molecular weight gradient (V^) is determined by 
the evolutionary state, but it can be also modified by differ- 
ent mixing processes taking place in the radiative interior. 
If the N profile changes because of a different V/i, then we 
could expect a signature in the avoided crossings frequencies. 

In each of the evolutionary states a - d we considered, 
the different mixing processes we analyse (classical and diffu- 
sive overshooting, microscopic and turbulent diffusion) pro- 
vide models with a central structure showing a different V/i 
and therefore different A'' profiles. In Fig. [10] we plot the 
hydrogen mass fraction and A'^ profiles in the central re- 
gion of 12 Boo A models in the evolutionary states 6 and c 
resulting from the calibrations with classical overshooting, 
with (short-dashed lines) and without microscopic diffusion 
(solid lines), and calibration with turbulent diffusion (long- 
dashed lines) . Though the evolutionary state is the same, the 
difference of chemical composition gradient (Fig. [lO] upper 
panels) is clearly refiected on the Brunt- Vaisala frequency 
profile (Fig. I fol lower panels), and therefore on the frequency 
of g-modes and on the properties of the avoided crossings 
(see Fig. Ill|l . Consequently, the determination of mixed- 
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Figure 10. Hydrogen profile (upper panel) and Brunt- Vaisala frequency {lower panel) in models computed witii overshooting (continuous 
line), with turbulent mixing (long-dashed line) and with diffusion of H and He (short-dashed lines) 



mode frequencies with a precision better than ~ 2/iHz for 
frequencies around 650 /xHz, would aUow to determine the 
properties of the chemical composition profile and, hence, 
of the physical mechanism responsible for the extra-mixing 
outside the core. The effect of different Vpi is even larger 
at lower frequencies: at 500 yitHz, where the expected mode 
amplitude is slightly larger than the half of the amplitude 
at fmax, the difference between overshooting and Dt models 
reaches 5 fjHz. 

We notice that such a valuable inference on the proper- 
ties of the chemical composition gradient is made possible, 
in the example presented here above, by a trade-off between 
two effects. On one hand, the difference in V^, which is 
due to a different mixing process, decreases as we consider 
more evolved models. On the other hand, the number of 
g modes entering the domain of solar-like oscillations, and 
therefore the number of avoided crossings, increases as the 
model evolves. This is the reason why, even though the dif- 
ferences between chemical composition profiles provided by 
classical overshooting and by Dt models increase as we con- 
sider less evolved models (see Fig. 1101 lower panel), the lack 
of £ = 1 avoided crossing predicted in the domain of solar- 
like oscillations for the scenarios c and d, makes difficult to 
infer information on the detailed properties of the /i gradi- 
ent. 

As mentioned above, models in the second gravitational 
contraction present £ — 2 avoided crossings in the solar-like 
spectrum domain. These 1 = 2 modes are however less sen- 
sitive than 1=1 modes to the A'^ profile, and the difference 
between overshooting and microscopic diffusion profiles is 



not large enough to be reflected in £ = 2 frequency modes. 
Models in the scenario c including turbulent diffusion show, 
however, very different N in the central region, in particu- 
lar, the model reaches the same location in the HR diagram 
with a smaller convective core, and this is reflected on a 
difference in the frequency of modes undergoing an avoided 
crossing, of the order of 100 /iHz with respect to those in 
the overshooting models. 

Concerning the main-sequence calibrated models (both 
for A and B components) we note that while in overshooting 
models 5vo2 as a function of v has a slope increasing with 
the value of the overshooting parameter, for turbulent dif- 
fusion models, 5uo2 shows only a weak //-dependence. That 
can be explained in terms of a different effect on the size of 
the convective core. While, in fact, in overshooting models 
the increase of aov implies a larger core, in turbulent dif- 
fusion ones the increase of Dt parameter decreases V/i at 
the border of the convective core, but leaves the core radius 
almost unaffected. Besides, the convective core for 12 Boo 
stellar masses is too small, and the difference between 5vo2 
slope for overshooting and Dt models not large enough to 
be detected in the expected solar-like frequencies. 

The effect of other diffusion processes such as the dif- 
fusive overshooting or the microscopic diffusion on the ^- 
proflle and on the frequencies of the avoided crossings are 
similar to those described for the turbulent mixing. The ef- 
fect of microscopic diffusion is, however, not limited to the 
core: it also modifies the structure of the envelope by a de- 
crease of the helium abundance in the outer convective zone 
and by the appearance of a chemical composition gradient 
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Figure 11. 1 = and £ = 1 large frequency separation in mod- 
els of which chemical gradient and Brunt- Vaisala are plotted in 
Fig. 1101 Model Ab with overshooting qqv = 0.15 (continuous 
line), and model with turbulent mixing (long-dashed line). The 
different chemical composition gradient (see Fig. 1101 affects the 
frequencies of £ = 1 modes undergoing avoided crossings (lower 
curves) but not the frequencies of radial modes (upper curves). 



below the convective envelope. The combinat i on of these ef- 
fects leaves, as noticed e.g. in iTheado et aP (|2005l ). a clear 
signature in the frequencies of acoustic modes. In order to 
fully exploit this seismic signature a high frequency reso- 
lution, that only long and uninterrupted observations can 
provide, is however needed. 



5 CONCLUSIONS 

We have presented a detailed modelling of the binary system 
12 Bootis by fitting the available observational constraints: 
effective temperatures, luminosities, chemical composition, 
and high precision masses of both components. As result 
of different 12 Boo calibrations we found a set of possible 
theoretical scenarios where the secondary is on the main 
sequence, whereas the primary has already left it. Its pre- 
cise evolutionary state, however, may vary from the SOC 
(models Ac) to a thick-shell-H-burning phase (models A a) 
by varying the overshooting parameter qqv from 0.23 to 
0.06. These values slightly decrease (by 0.02) if microscopic 
diffusion is included in stellar modelling. Other central mix- 
ing processes, such as diffusive overshooting (described by 
the P parameter), and a turbulent mixing with a paramet- 
ric turbulent diffusion coefficient (Dt) lead to the similar 
results. For each transport processes we found a possible 
range of values of the parameters aov , P and Dt that place 
12 Boo B in the MS and 12 Boo A between the second 
gravitational contraction and the sub-giant phase. The only 
way to place both components in the MS (models Ad) is to 
assume a different efficiency of the mixing process in both 
components, for instance : aov, a = 0.37, aov,B = 0.15, or 
-Dt.a = 330 cm^ s"\ L>t,b = 100 cm^ s"^ With the avail- 



able observational constraints we are not able, however, to 
discriminate among the different scenarios and models pro- 
posed: it is cleat that additional and independent observa- 
tional constraints are needed. 

In Sec.|4]we have shown that the detection of the theo- 
retically predicted solar-like oscillations in 12 Bootis would 
provide a powerful test to discriminate between different 
scenarios. Indeed, the detection of avoided crossings in the 
primary oscillation spectrum will allow a robust inference 
on the evolutionary state of the star and, therefore, on the 
amount of overshooting (or extra-mixing) from the core. The 
presence or absence of ^ = 1 mixed modes will allow to dis- 
criminate between sub-giant and MS models, and the ab- 
sence of both 1 = 1 and 1 = 2 will allow to place 12 Boo A 
in the MS, before the SOC. This possibility is especially ex- 
citing since it would imply that, even for stars with masses 
so close, the central mixing has might have been working 
with a significantly different efficiency, opening the door to 
important questions concerning transport processes. 

As shown in Sec. 14.21 for a model in a given evolu- 
tionary state, a different chemical composition profile in the 
core can be due to a different physical process responsi- 
ble for the extra-mixing. This affects the behaviour of the 
Brunt- Vaisala frequency and hence the frequency of modes 
of mixed p- and g- character. Consequently, information on 
the mixing process taking place in the stellar center could 
be inferred from the mixed mode frequencies. The precision 
required to reach this goal is, however, evolutionary state de- 
pendent. While for a thin-shell-H-burning phase, a precision 
of 2 /iHz should be enough to distinguish between the sharp 
chemical profile provided by overshooting and the smother 
one from a diffusive process, in a slightly more evolved phase 
(thick-shell-H-burning phase) the effect would be already 
erased. On the other hand, as stellar age decreases, the num- 
ber of mixed modes also decrease. As consequence, for mod- 
els in the second gravitational contraction, information on 
the mixing process can be derived from 1 = 2 mixed modes 
only if the generated chemical profiles are quite different. For 
MS models, where no avoided crossings are predicted, the 
effect of the central mixing on the small separation is 
not large enough to distinguish between different transport 
processes. 

Finally, the detection of solar-like oscillations also in 
the B component, will reduce the uncertainties on the HR 
location of 12 Boo B, and hence will significantly improve 
the observational constraints of the system. Nevertheless, 
as all the models of component B are in the MS, no useful 
information concerning the age of the system will be add by 
the 12 Boo B frequencies. 

In conclusion, we have shown that thanks to the precise 
knowledge of the masses of 12 Bootis, and the additional 
constraints due to the binarity, the search and detection of 
solar-like oscillations in 12 Boo, hopefully by a spectroscopic 
multi-site campaign of observations, would mean an impor- 
tant step in the understanding of stellar evolution in the 
mass domain where the the convective core appears. 
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